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RADIATION GEOMETRY FACTOR BETWEEN THE EARTH AND A SATELLITE 


SUMMARY 


In solving thermal problems of satellites, the geometry factors must be known. 

A brief discussion of geometry factors, an analytical expression, and computer procedure 
is given. 

Methods for computing these geometry factors from basic relations require 
lengthy computer programs with many inputs. A method was devised in the Space Thermo- 
dynamics Branch of the Research Projects Laboratory of MSFC for obtaining these 
factors for space thermal calculations from a subroutine similar to the common trigono- 
metric subroutines. It requires only 600 octal locations. This is less than 10 percent 
of the total memory of the IBM 7094 Mod. II used by the Computation Laboratory. Max- 
imum computer time required per factor obtained is . 0015 seconds on the IBM 7094 Mod. 

II. This subroutine has been used for over two years without any apparent problems. 


INTRODUCTION 


The geometry factor, F, between two bodies of various geometrical configurations 
must be known in order to solve most heat problems encountered in space. A brief de- 
scription of geometry factors follows. 

The Stefan-Boltzmann radiation law (fourth power) relates the amount of radiative 
flux emitted from a black body to its absolute temperature. A non-black body radiates a 
fraction of the amount of energy flux radiated by a black body at the same temperature. 
This fraction is known as the total hemispherical emittance (or emissivity) of the non- 
black body at the given temperature. 

Assuming opaque bodies and Lambert's cosine law, the angular distribution of the 
flux radiated from an elemental area on a body, (a) , can be calculated. If another body, 
(b) , is placed in the vicinity of (a) , each elemental area on (b) has incident upon it some 



fraction of the flux from each elemental area on ( a) . It follows that the amount of incident 
flux absorbed on (b) can be calculated by a double integration. This amount is a fraction 
of the total flux emitted, and this fraction is called the geometry factor, F. 

Reflected parallel radiation from an outside source, or albedo, also enters into 
orbital problems. Assuming diffuse reflection, an albedo-F can be obtained by a similar 
double integration. Except for the simplest geometrical configurations, the integrations 
for F and albedo-F must be carried out by numerical techniques. 

A fact which is often overlooked is that, by using the reciprocity theorem, the 
F , -r, (the geometry factor from a body, A, to another body, B) can be calculated if the 

A *15 

F_ . (the geometry factor from the body, B, to the body, A) is known. Obviously, the 
B— 

simplest route in computation should always be chosen. For example, the F is much 
more easily calculated between the earth and a smaller object, such as a satellite, by 
thinking in terms of the satellite radiating to the earth instead of vice versa,, although in 
practical problems the F is used to compute earth radiation to a satellite. 

Convair Corporation performed the necessary integrations, previously mentioned, 
for several geometrical configurations. Their results were published in graphical form 
[ 1 , 2 ]. 


Presently, the Space Thermodynamics Branch of Research Projects Division is 
using these results for two of these geometrical configurations for both albedo and the so- 
called IR (earth's infrared) , which is the earth's Stefan- Boltzmann radiation. The two 
configurations are: ( 1) a sphere (earth) and a flat plate (on a satellite) , and (2) a sphere 
( earth) and a cylinder ( on a satellite) . 

Since most thermal problems of orbiters require many tedious arithmetic opera- 
tions, the IBM 7090 computer in the Computation Division has been utilized. An analyti- 
cal expression for F as a function of the known geometrical parameters (such as attitude, 
etc.) is imperative for computing purposes. Calculating the factors from basic assump- 
tions by numerical integration is too lengthy to use for a subroutine as visualized by R- 
RP-T. Since polynomials are very adaptive to computers in that they require a minimum 
of memory storage, minimum computing time, and are easily programmed, an empirical 
polynomial was generated with a "sum of the least squares polynomial curve fitting" 
routine at the Computation Laboratory. 

PROCEDURE 


A. F FOR EARTH'S IR TO CYLINDER 

The following parameters are geometrically defined in figure 1 for IR be- 
tween a planet and a cylinder: 

2 



r = altitude (km) 
y = attitude angle ( deg) 


Figure 2 shows the F versus altitude for a cylinder as plotted in [ 2] for several 
representative attitudes. Figure 3 shows the 'same curves as figure 2 on Cartesian 
coordinates over the range 182 km :=r^ 3500 km for y = 0° , 20° , 40° , 60°, and 90°. 

Letting F be a function of h for each curve in figure 3, a least square curve fit 
equation was obtained of degrees 2 through 9. The 5th degree equations were determined 
to be the smallest degree acceptable. Thus, for each y value (0®, 20°, 40° , 60° , 90°) , 
a 5th degree polynomial equation was obtained: 
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Other least square curve fit processes were used to obtain the as a function of 
y. Of these, the 4th degree polynomial equations were best: 


B.. 

1J 


Z 


k=0 



9 


with the A., given in table II. Hence, 
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for a cylinder 0° < y^_ 90°, since the cylinder is symmetrical about its axis. The above 
equations are not valid outside the limits 0° ^ y ^ 90° and 182 km ^ r ^ 3500 km. 

B. F FOR EARTH'S IR TO FLAT PLATE 


Similar figures (4, 5 & 6) and tables (III & IV) are shown for a flat plate 
radiating on one side. The resulting equation for the F of a flat plate is: 


F 



Z 


k=0 



k 

y 


i 


r 


where 0° s y < 110° and 182 km ^ r ^ 3500 km. 


Since y can range from 0° to 180° for a flat plate, and the case 110° to 180° is not 
covered in the polynomial, most of the F's in this range are zero for all practical pur- 
poses. The non-zero F's are obtained by a linear interpolation included in the subroutine 
for F. The interpolation routine is not as adequate as desired, due to the fact that F is 
not a linear function of y or h and F-*0 at different h values for a given y. This fact makes 
it necessary to interpolate differently over specified increments of altitude. This causes 
a discontinuity in interpolative values of F at the value of h where any two interpretative 
groups are separated. A discontinuity (less than ± 3% F max ) in computation occurs upon 
crossing one of the boundaries. It should be noted that this occurs when 0 < F < . 2 . 


C. ALBEDO-F FOR BOTH CYLINDER AND FLAT PLATE 


Figures 7 and 8 show a graphical definition of the parameters necessary to 
compute the albedo-F. 

The albedo-F is not only a function of parameters y and h, but of 0 and <f>_ also. 

O L/ 

These angles did not enter into the IR case because of symmetry. Figures 9 through 18 
(taken from [2] ) are typical examples of graphs for albedo-F. 

Intuitively, the dependence of the albedo-F upon 0 S is thought to be approximately 
proportional to cos 9 S and that F = albedo-F when 0 g = 0. Also, albedo-F is a week 
function of < p c . These assumptions are illustrated in figures 17 and 18 where the solid line 
curves are typical integrated curves and the dotted line curves are the corresponding 
curves obtained by using the above assumption. Hence, 

albedo-F = F cos 9 , if cos 9 > 0 

s’ 

and albedo-F = 0, if cos 0 < 0 . 

This subroutine is coded to calculate F for the flat plate or the cylinder. Only F 

is in the output; if the albedo-F is desired, one must multiply by cos 0 . 

s 
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CONCLUSION 


As far as the IR is concerned, there is one basic type of error (neglecting the 
interpolation routine). Assuming the data from the Convair reports to be valid, this is 
the F obtained for values of y not used in the curve fitting process. Since the By and Ajj 
are continuous, the error is estimated to be within ±2% of F max . Over an entire orbit, 
it is believed that these errors will usually have a cancelling effect. For the albedo-F the 
same errors occur as in IR, plus an error, brought about by using the assumption albedo- 
F=F,cos 0 g , when cos 0 g > 0 and albedo-F = 0, when cos 9 S 0, which is estimated to be 
less than ± 5% of F 

max 

Methods for computing these geometry factors from basic relations require 
lengthy computer programs with many inputs. A method was devised in R-RP-T for 
obtaining these factors for space thermal calculations from a subroutine similar to the 
common trigonometric subroutines. It requires only 600 octal locations. This is less 
than 10 percent of the total memory of the IBM 7094 Mod. II used by the Computation 
Laboratory. Maximum computer time required per factor obtained is . 0015 seconds on 
the IBM 7094 Mod. II. This subroutine has been used for over two years without any 
apparent problems. 

In making typical computations of the radiation, it must be remembered that the 
factors are defined such that one uses the actual area of one side of the plate and the 
projected area (length x diameter) of the cylinder. 
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FIGURE 2. ■ GEOMETRIC FACTOR FOR EARTH THERMAL RADIATION INCIDENT TO A 
CYLINDER VS. ALTITUDE, WITH ATTITUDE ANGLE AS A PARAMETER 
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FIGURE 5. GEOMETRIC FACTOR FOR EARTH THERMAL RADIATION INCIDENT TO A 
FLAT PLATE VS. ALTITUDE WITH ALTITUDE ANGLE AS A PARAMETER 
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FIGURE 6. GEOMETRY FACTOR VS. ALTITUDE FOR IR TO A FLAT PLATE 
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FIGURE 9. GEOMETRIC FACTOR FOR ALBEDO TO A CYLINDER VS. ALTITUDE WITH 
ZENITH DISTANCE BETWEEN CYLINDER AND SUN AS A PARAMETER 
y = 0°, <£ c = 0° - 180° 
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GEOMETRIC FACTOR 
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FIGURE 10. GEOMETRIC FACTOR FOR ALBEDO TO A CYLINDER VS. ALTITUDE WITH 
ZENITH DISTANCE BETWEEN CYLINDER AND SUN AS PARAMETER 
7 = 60° , <$> c = 0° 
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GEOMETRIC FACTOR 



FIGURE ii„ GEOMETRIC FACTOR FOR ALBEDO TO A CYLINDER VS. ALTITUDE WITH 
ZENITH DISTANCE BETWEEN CYLINDER AND SUN AS A PARAMETER 
y = 60 ° , <p c = 90 ° 
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GEOMETRIC FACTOR 
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FIGURE 12. GEOMETRIC FACTOR FOR ALBEDO TO A CYLINDER VS. ALTITUDE WITH 
ZENITH DISTANCE BETWEEN CYLINDER AND SUN AS A PARAMETER 
y = 60° , (p c = 180° 
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FIGURE 13. GEOMETRIC FACTOR FOR ALBEDO TO ONE SIDE OF A FLAT PLATE VS. 

ALTITUDE WITH ZENITH DISTANCE BETWEEN FLAT PLATE AND SUN 
AS PARAMETER y = 0° , <£ c = 0° - 180° 


18 



GEOMETRIC FACTOR 



IOO 500 1000 5000 

ALTITUDE-nautical miles 


FIGURE 14. GEOMETRIC FACTOR FOR ALBEDO TO ONE SIDE OF A FLAT PLATE VS. 

ALTITUDE WITH ZENITH DISTANCE BETWEEN FLAT PLATE AND SUN 
AS PARAMETER y = 30°, <p c = 0° 
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FIGURE 15. GEOMETRIC FACTOR FOR ALBEDO TO ONE SIDE OF A FLAT PLATE VS. 

ALTITUDE WITH ZENITH DISTANCE BETWEEN FLAT PLATE AND SUN 
AS PARAMETER y = 30° , <£ p = 90° 
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GEOMETRIC FACTOR 



FIGURE 16. GEOMETRIC FACTOR FOR ALBEDO TO ONE SIDE OF A FLAT PLATE VS. 

ALTITUDE WITH ZENITH DISTANCE BETWEEN FLAT PLATE AND SUN 
AS PARAMETER y = 30° , <f> c = 180° 
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FIGURE 18. GEOMETRY FACTOR VS. ALTITUDE FOR ALBEDO TO A FLAT PLATE 
(7 = 30 °) 
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TABLE I 
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( 7 = 60=) 
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+. 13303x1 0 1 
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+. 14086x1 0 1 

t— <- 

II 

-.U546xi0" 2 

12129x1 0“ 2 
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95335x1 0 " 3 
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i = 2 
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+.58425x10 6 - 
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+.739llxl0 _13 

+.73500xl0‘ 13 

+. 46419xl0~ 13 

+.3 0839x1 0" 13 
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A.. FOR CYLINDER 
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A.. FOR A FLAT PLATE 
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-. 12748xl0 -5 

+, 14561x10 
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